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The cosmological model consisting of an electromagnetic Born-Infeld (BI) field coupled to a 
Robertson- Walker geometry is tested with the standard probes of SNIa, GRBs and direct Hub- 
ble parameter. The analysis shows that the inclusion of the nonlinear electromagnetic component 
does not contribute in a significative way to the observed expansion. The BI electromagnetic matter 
is considered with an abundance of Qbi, that our best fit leads to Q,bi = 0.037 when tested with 
SNIa and the Hubble parameter data (0.1 < z < 1.75); while when tested with GRBs the result 
is of Q.bi = 0.304, which may indicate that this electrodynamics was important at epochs close to 
the appearance of large structure (z « 7), although this late result has not as much reliability as 
that corresponding to the first two probes, since we know that the dispersion in GRBs data is still 
considerable. In view of these results we can rule out the electromagnetic Born-Infeld matter as 
the origin of the present accelerated expansion, this conclusion concerns exclusively the Born-Infeld 
theory. 



I. INTRODUCTION 

In the quest to find out the dark energy nature, there have been many reasonable proposals, other than the 
cosmological constant, some of them modeling dark energy very succesfully. Inspired in some inflationary models, 
nonlinear electromagnetic matter has been invoked lately to mimic the dark energy behaviour. 

Nonlinear electrodynamics (NLED) has been addressed to produce inflation in the early universe [UH], and it has 
also been associated to the magnetic seeds that generated the large scale correlated magnetic fields observed nowadays 

[3H51- 

According to Einstein's equations and assuming a Robertson- Walker (RW) geometry, the accelerated expansion is 
attributed to a kind of repulsive gravity effect. Such expansion is only possible if the dominant component of the 
universe has effectively a negative pressure that counteracts the attractive effect of ordinary matter, with the energy 
density p and pressure p being such that p + 3p < (so the strong energy condition (SEC) is violated). 

It has been shown that nonlinear electrodynamics can violate the strong energy condition, consequently producing 
negative pressures that in turn accelerate the expansion. In fact, several models have been proposed that consider 
nonlinear electromagnetic matter to produce accelerated expansion [SHH] ; and so these models open the possibility 
that dark energy is rooted in properties of ponderable fields and matter. 

The violation of the SEC can also be associated to a nonvanishing trace in the energy momentum tensor, T — 
T£ = g^Tpw With the apropriate sign, the contribution of pressure T can provide the anti-gravity effect associated 
with dark energy, in this sense the phenomenologically nonvanishing NLED energy-momentum trace is worth to be 
considered [5]. 

The similarity of T to the cosmological constant was already noted in [10]: focusing on the trace anomaly of 
quantum chromodynamics, the issue of whether nonperturbative effects of self-interacting quantum fields in curved 
spacetimes may yield a significant contribution was addressed in jTUj , where a preliminary estimate of the expected 
order of magnitude indicated the potential relevance of this effect. 

A similar argument is that General Relativity can be viewed as a low energy quantum effective field theory of 
gravity, provided that the Einstein-Hilbert classical action is endowed with the additional terms required by the trace 
anomaly [llj . 

There is another approach that considers a nonlinear behaviour in the propagation of light, similar to light traveling 
in non vacuum spacetime (cosmic substratum) |12j : this approach is supported by the form of the nonlinear electro- 
magnetic Born-Infeld equations that are formally identical to Maxwell's for a substance with a dielectric permeability 
and magnetic susceptibility being certain functions of the field strength but without a spatial distribution of charge 
and current [T3]. The Born-Infeld theory, a classical theory designed to regulate point-particle-induced divergences, 
provides a more complete theory of electromagnetism, that includes nonlinear quantum electrodynamic (QED) effects 
as vacuum polarization. Moreover, it is now also considered as an effective theory arising from string theory |14j . The 
breaking of scale invariance in the Born-Infeld (BI) action has the immediate consequence that the resulting equation 
of state allows for negative pressure. 
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Producing accelerated expansion is not enough to consider such proposals as an explanation for dark energy; 
the corresponding models should produce certain precise amount of acceleration and at a precise time. Therefore, 
confrontation with observations is mandatory to determine to what extent these models are able to reproduce the 
observable data or at least the lately observed accelerated expansion. In this paper we investigate the observational 
constraints on the homogeneous and isotropic cosmology coupled to an isotropic electromagnetic field governed by a 
BI lagrangian. 

The Robertson- Walker (RW) spacetime describe a homogeneous and isotropic geometry, which is suggested by 
cosmological observations, however coupling the RW spacetime to an electromagnetic field originates a main com- 
plication: electromagnetic field single out spatial directions breaking the isotropy of RW spacetime. One way to 
overcome this difficulty is to adopt a spatial average in the electromagnetic field, [MS]) [T51 116j . 

Alternatively we can consider a vector triplet compatible with space homogeneity and isotropy of RW [TTl [TB] . 
This is a set of three equal length vectors that point in three mutually orthogonal spatial directions. While the triad 
guarantees the isotropy of the background, it does not automatically imply the isotropy of its perturbations, as it 
could be necessary to model some observed anomalies in the CMB radiation. In fact, the cosmic triad can be realized 
with a classical SU(2) vector field configuration [T71 [T§1 |2"U] . 

In Sec. II we couple the BI electromagnetic field to FRW. In Sec. Ill we introduce the SNIa, Hubble parameter 
and GRBs probes as well as the statistical method. In Sec IV we present the constraints obtained as well as the cases 
in the Born-Infeld literature that are included in our analysis; in the last section we draw some conclusions. 



II. COUPLING BI TO FRW 

We assume the homogeneous, isotropic, spatially flat geometry that observations led to conclude 

ds 2 = dt 2 - a 2 (t){dr 2 + r 2 d9 2 + r 2 sin 2 9dcj) 2 ), (1) 

where a(t) represents the scale factor of the Universe. 

In this context the dynamics of the Universe can be derived from the Friedmann equations, once the characteristics 
of the fluid that models the cosmic substratum have been established. In the present case, besides the dark matter 
and dark energy, we include a nonlinear electromagnetic fluid component characterized by a lagrangian L(F, G) 
that depends nonlinearly on the invariants of the electromagnetic field, F = F^F^" and G = F^F^" , and F^ — 
( y /—g/2)e fJiVpcr F p ' 7 is the dual to the electromagnetic field tensor F„„. 

To plug in an electromagnetic field to the isotropic RW geometry, an average process is needed to isotropize 
the preferred-directions of the electromagnetic field; afterward we can consider the electromagnetic field as a fluid 
characterized by isotropic pressure p and energy density p. 

When working with a U(l) field, several authors adopt the spatial average proposed by Tolman and Ehrenfest 
(1933) [T5], that consists in considering the thermal equilibrium of a general static gravitational field which could 
correspond to a system containing solid as well as fluid parts. This average can be posed as follows [16]: the spacetime 
is filled with electromagnetic radiation, the one with cosmological interest is that of the CMB. It can be considered as 
a stochastic background of short wavelength radiation (compared with curvature) that satisfies the spatial average for 
volumes that are large compared to the wavelength but small compared with the curvature (R » volume >> Xcmb)- 
Other authors 6 consider that the isotropization of the electromagnetic field occurs at the quantum level and arises 
from the condensates that appear due to vacuum fluctuations. 

Alternatively some authors assume the vector field as part of a cosmic triad [171 118j i.e. a set of three identical 
vectors pointing in mutually orthogonal spatial directions. 

The cosmic triad can appear naturally from a non-Abelian SU(2) gauge theory, and some authors have adopted 
this realization of the triad [T7], [15] j f° r instance in [5DJ the non-Abelian Born-Infeld field in RW space was studied 
and the field equations were decoupled, deriving analytical expressions for the scale factor and the electromagnetic 
scalar function that determines the cosmic triad. 

Taking any of the above procedures, we end up with an isotropic energy-momentum tensor, with energy density 
p = Tq and pressure p = — T?/3, i = 1,2, 3, that for a nonlinear lagrangian L(F, G) are given by 

T^v = -4L F F™F av + (GL G - L)g^ = (p + p)u^u v -pg^ v , 
p = —L + GLq — AE 2 Lp, 

p = L-GL G + ^(E 2 -2B 2 )L F , (2) 

where Lx = dL/dX. 



3 



We shall consider the Born-Infeld action, coupled to gravity: 

S = J V^gd 4 x [i? + 4Gft 2 (Vx - 1) 



X — 1 + _L F F^ — (F F^") 2 

_ F G^_ 

~ 2/32 16/?4 ' ( 3 ) 

where R is the scalar curvature and ft is the BI critical field strength. Our convention for the Maxwell lagrangian is 
^Max = —F/4. 

Coupling to a RW geometry, from the Friedmann equations we get 

a \nG , „ , 

- = ^(p + 3p), 

a 3 

(4) 

where the overdot means differentiation with respect to the cosmic time t, and we have put c = 1. In addition the 
expression for the Born-Infeld energy density p = ps/, is given by 



Pbi = ft 2 (/X - l) . (5) 



From the previous equation, along with the energy conservation law, = p + 3H(p + p), the scaling between the 
electromagnetic fields and the scale factor a(t) can be obtained. To this end we take in particular X = 1 + F/2ft 2 , 
i.e. G — 0. Then we find that F = 2{B 2 — E 2 ) = (const)a -4 , consequently, electric and magnetic fields scale as 
E,B ~ a~ 2 . This result does not depend on the particular analytic dependence of L(F). 

The election of the constant that relates the proportionality determines the energy density sign, 



Pbi 



= ft 2 yjl - (const)a- 4 - 1 . (6) 



We can mention three possible cases: 

(i) If < (1 — (const)a -4 ) < 1 and (const) > 0, then psi < and matter is phantom-like, without negative 
kinetic energy. 

(ii) If < (1 — (const)a~ 4 ) and < const < 1, then psi may be greater or less than zero. 

(iii) If (1 — (const)a~ 4 ) > 1 and const < 0, then p BI > 0. 
Different authors chose different constants (see section IV. D). 

Considering the previous equations and that in a RW cosmology, a = (1 + we can write the expression for the 
Hubble parameter as a function of z and then proceed to determine to what extent the Born-Infeld electromagnetic 
field contributes to the observed accelerated expansion of the nowadays observed Universe. In addition to the BI 
electromagnetic field contribution, we have included a cosmological constant and matter (dust) to ensure cosmic 
acceleration. The observational tests have the aim to determine if there is a significant contribution from the Born- 
Infeld electromagnetic field, that shall be reflected in a lower contribution of the cosmological constant to produce 
accelerated expansion. 

From the Friedmann equation, the Hubble parameter in terms of the fractional energy densities and the redshift is 



H 2 (z) 



= n, 



«(i + zf + n A + n BI (yi + 3(1 + z) 4 - i) , (7) 



where flsi — Zg/SHQ with g — AirGfi 2 is the coupling constant of the BI field. We have fixed the constant in Eq. 
as const= —3 [20 . In this way, the matter components satisfy the following normalization condition 

n m + n A + n BI = 1. (8) 

Notice that the BI field scales the Hubble parameter, roughly speaking, as (1 + z) 2 , so from the very beginning one 
can predict it will not be very important at present time. The corresponding density of the BI component falls off 
more slowly than radiation, H 2 ~ (1 + z) 4 , and presureless matter, H 2 ~ (1 + z) 3 , probably was significative in far 
epochs as reionization. Moreover, from our analysis we shall obtain an estimation for the present Born-Infeld field ft. 

In the following section we briefly describe the used statistics to compare the models with the most recent SNIa, 
direct Hubble parameter and GRBs data. 
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III. OBSERVATIONAL DATA AND STATISTICAL METHOD 

We derived the constraints to the model's free parameters using three different observational data sets: 

• the Union2.f supernovae data set, the updated of Union2 reported in [21] . 

• the direct Hubble parameter data given in |22j . 

• the Gamma-ray Bursts luminosity distances calibrated in (23) . 

A. Supernovae: N. Suzuki et al. 2011 data set 

Recently, the Supernova Cosmology Project (SCP) collaboration released the updated Union2.1 compilation which 
consists of 580 SNIa [21]. The Union2.1 compilation is the largest published and spectroscopically confirmed SNIa 
sample to date. The data points of the 580 Union2.1 SNIa compiled in [21] are given in terms of the distance modulus 
/J-obs(zj). On the other hand, the theoretical distance modulus is defined as 

H(zj) = 5log 10 [d L (zj,6i)] +fi , (9) 

where di,{ZjA) is the Hubble free luminosity distance 

c (i + z) r* dz> 

dL{zA) = -^Jo Way (10) 

and 6i denotes the model parameters. 
The x 2 function for the SNela data is 

2/,, a\ (jJ^jjJMhQd ~ fJ-obs(Zj)) 2 . . 

x M (Mo»f<j = 2^ jh^a ' " x > 

where the a 2 ^ corresponds to the 68.3% error of distance modulus for each supernova. The parameter /xo is a nuisance 
parameter that encodes the Hubble parameter and the absolute magnitude M, and has to be marginalized over. Here, 
we follow the method suggested in [24, 25 which consists in maximizing the likelihood by minimizing x 2 with respect 
to /Wo- Then one can rewrite the \ 2 as 

Xsn(8) =c%- 2c 2 /io + c 3 Mo, (12) 

where 

C1 ~U ' (13) 



\- (fifa; fj, = 0, 6i) - (j, obs (zj)) 



02 -2^ ' (14) 



580 

C3 = 



(15) 



The minimization over /i gives fig — C2/C3. So the \ function takes the form 

XSN(8i) =c x -^. (16) 
C3 

Since Xsn = Xsn(^o = 0: $i) (up to a constant), we can instead minimize Xsn which is independent of /zq. 
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B. Hubble parameter: M. Moresco et al. 2012 data set 



We use the compilation of Hubble parameter measurements estimated with the differential evolution of passively 
evolving early-type galaxies as cosmic chronometers, in the redshift range < z < 1.75 recently updated in |22j but 
first reported in [26] . 

The root idea supporting this approach is the measurement of the differential age evolution of these chronometers 
as a function of redshift, which provides a direct estimate of the Hubble parameter H(z) = — 1/(1 + z)dz/dt ~ 
— 1/(1 + z)Az/ At. The main strength of this approach is the confidence on the measurement of a differential quantity, 
Az/At, which provides many advantages in minimizing many common issues and systematic effects. In addition, 
compared with other techniques, this approach provides a direct measurement of the Hubble parameter, and not of 
its integral, in contrast to SNela or angular /angle- averaged BAO. 

The best fit values of the model parameters from observational Hubble data are determined by minimizing the 
quantity 

Xh( H 0,Vi) - „2 71 \ ' i 17 ) 

J = 1 a H\ z j) 

where are the measurement variances, H Q will be fixed from |27) . while the vector of model parameters, 9i, will 
be Oi = (f2 m ,f2 S j) after using Eq. (|]) in Eq. @- 



C. GRBs: Wang 2008 data set 



If one is interested in knowing how dark energy behaves, one way of improving our understanding is through the 
study of the evolution of the Universe at redshifts higher than those probed by SNe la, then one has to resort to 
visible objects at large distances, if any. 

So far, GRBs are the most luminous astrophysical events observable today and can be seen up to z ~ 8 28J. 
However, the degree to which GRBs can be used as standard candles is not yet fully understood, and to extract 
cosmological information, it is necessary to calibrate them for each cosmological model tested. 

Among the methods to calibrate GRBs in a cosmology-independent way, there are proposals that use SNe la 
measurements to calibrate GRBs externally, so they can be used to constrain dark energy models 29 31 . Indeed, an 
update has been made recently of the sample of GRBs now consisting of 109 objects [32] ■ In [52J[33] new calibrations 
have been developed in order to use GRBs in cosmological tasks. 

Here it is worth noting that cosmological constraints from GRBs are sensitive to how these are calibrated. Calibrat- 
ing GRBs using Type la supernovae gives tighter constraints than calibrating GRBs internally [33] ■ In the present 
work we have chosen to work with the calibrated sample in |23j . The advantage of this sample is that it is internally 
calibrated, and then it may be straightforwardly combined with other data when deriving cosmological constraints, 
which is actually our aim. However, we have to point out that this method also requires a cosmological model as an 
input and thus it is not completely cosmology independent. 

The x 2 function for the GRBs data is 



Xgrb = [Ar p (zi)] • (Cov G Jj B ) . • [Ar p ( 



PV'J ) 



(18) 



Ar p (zi) 



-data i 



rp"- J (zi) -fp(zj), 



where r p (z) is defined by 



r p {zi) 



r p( z ) 
r„(0.17)' 



(l + z)^H 



r(z), 



and r(z) = cIl(z)/(1 + z) is the comoving distance at z. 
The covariance matrix is given by 



(Cov GRB ) i:j = a(r p (zi))a(r p (zj)) (Cov G rb) 1 



(19) 



(20) 



where Govgrb is the normalized covariance matrix given in Table 3 from [23] . and 

a{f p { Zi )) = a (r p (z, ; )) + , if r p (z) > r p (z) data ; 
a(r p (zi)) = cr [r p { Zl )y , if r p {z) < r p (z) data , 



(21) 



() 



1.0 



0.6 



0.4 



0.2 



0.0 



0.1 




0.4 



0.14 
0.12 
0.10 
0.08 
0.06 
0.04 
0.02 



\ 



o.oo' 

0.10 0.15 0.20 0.25 0.30 0.35 0.40 



(b) 



Figure 1. Confidence regions in the (Q m , Qbi) plane for the Born-Infeld cosmology from the different observational data sets, 
(a): green ones correspond to GRBs, blue to Hubble parameter data, orange to SNela. (b): blue to the combination of SNefa 
and Hubble parameter data and green to the combination of all observational data. To construct these confidence regions we 
have considered a prior on Sl m = 0.266 ± 0.029 from WMAP-7 years. 



where a (r^(zi)) + and a {f^{zi))~~ are the 68% C.L. errors given in Table 2 from |23j . 

In order to constrain the model parameters, we use a Markov Chain Monte Carlo (MCMC) code to maximize the 
likelihood function C(9i) oc exp[— x 2 (0i)/2] where 9i is the set of model parameters and the expression for x 2 {@i) 
depends on the dataset used. 

The MCMC methods (completely described in [51|-f3l)] and references therein) are well-established techniques for 
constraining parameters using observational data. To test their convergence we follow the method developed and fully 
described in [37]. 

On the other hand, it has been argued that the presence of a prior on f2 m leads to uniformity in its best fit 
value, thus minimizing its influence in the dark energy constraints |38j . Here, we choose a Gaussian prior on Q m — 
0.266 ± 0.029 from WMAP-7 years [3J5] and also for running our MCMCs we adopt the physical controls < O m < 1 
and < Q B i < 1. 



IV. RESULTS 



A. Cosmological Constraints 

Our constraints on model parameters are given in Table [i] and illustrated in Figure [T] The results for il m are 
entirely consistent with the amount of matter predicted by the standard model ACDM, so the discussion will focus 
on the values obtained for the nonlinear electromagnetic matter component. From Table [T] we can see that most of 
the constraints for model parameters, from SNela, Hubble parameter and combinations of them, suggest that the 
contribution of Ubi is weak, however the best fit value for Ubi obtained from the observational data of GRBs seems 
to be more significant. Besides, from Figure [TJa), we see that restrictions from Born-Infeld cosmology are stronger 
if we use SNela or H(z) in contrast to GRBs, and also from Figure [ljb), we see that the combination of SNela and 
H (z) data sets yields solid and strong constraints, but from the combination of all observational sets the confidence 
region seems to be slightly less tighter. 



B. Equation of state parameter 

The required gravitational properties of dark energy needed to induce the accelerated expansion are well described 
by its equation of state w(z) = p/p, which enters into the Friedmann equations (second in Eqs. Q), implying 
that a negative pressure (w < —1/3) is necessary in order to induce accelerated expansion. So, the parameter w(z) 
determines not only the gravitational properties of dark energy but also its evolution, which can easily be obtained 
from the energy momentum conservation. So, from a general expression of the Friedmann equation, Eq. Q, 
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Table I. Summary of the best estimates of model parameters. The errors are at 68.3% confidence level. See Figure [T] 



w(z — 0) dw/dz\ z=0 q z t 
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Table II. Summary of the best values for the equation of state parameter w(z — 0), the deceleration parameter qo and the 
transition redshift (deceleration/acceleration) Zt- The errors are at 68.3% confidence level. 
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Vl x exp 3 
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the knowledge of Sl m , z and H(z) suffice to determine w(z) as follows: 

2d\nH 

w(z) = — 



1 



3 dz 
Ho 
H 



(l + z)-l 



r» m (i + z)3 



Rewritting Eq. (231 in terms of the parameters model, Qbi and f2 m , we obtained 

2fl BI (l + z) 4 



w{z) 



-1 + 



v/l + 3(l + z) 4 1 - 2Q BI -Q m + S7 V 1 + H 1 + Z Y 



(22) 



(23) 



(24) 



Note that we have considered that Q\ and Qbi constitute an unified entity. 

In Table [IT] we summarize our results and in Figure [2] we present the evolution of w(z) using the best fits from 
each observational set. In general our results point to a quintessence behaviour with w(z) > —1 and also that 
dw/dz\ z= o > for the current time. As for the GRBs probe, the tendency of w(z) seems to be less negative. 



C. Deceleration parameter 



Regarding the so-called kinematic approach, in which the deceleration parameter q is parameterized as a function 
of the redshift z, it is straightforward to obtain q as function of the free parameters of the model using the equation 
of state parameter, Eq. (23), and the Hubble parameter H(z), Eq. Q, as follows: 



Q(z) 



1 - 



Qm (1 + z)" 
H 2 



w(z) 



(25) 




(a) SNela (b) Hubble (c) SNela+Hubble 




(d) GRBs (e) SNela+Hubble+GRBs 



Figure 2. Equation of state parameter w(z) using the best fit for Q m and Q.bi from each observational data set. The central 
line represents the best fit and the shaded contour represents the lcr confidence level around the best fit. 
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The different values for the evolution of the expansion rate are presented in Table [T] and also the results of the com- 
putation q(z t ) = to obtain the transition redshift, z t . All these quantities are given with uncertainties corresponding 
to the 68.3% interval of confidence. The evolution of q with z can be seen in Figure [3] 
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D. Born-Infeld cosmologies 



Our analysis includes and applies to, at least, the following works on Born-Infeld cosmologies: 

1) The one presented by Novello et al. [8] for Born-Infeld magnetic universe, with U(l) electromagnet ism and 
performing the Tolman spatial average; the energy density given by Eq. |5]), with X = 1 + F/(2(3 2 ), i.e. with 
vanishing electric field and therefore G = E ■ B = (cf. Eqs (B.2) and (B.3) in [§]). This scheme is very similar to 
that presented in [16] with the identification (3 2 — > (26 2 ) -1 . 

2) Dyadichev et al. [5D] considered a cosmic triad realized by a SU(2) field parametrized in terms of a function of 
time h(t) as 



E? = h{t)5t, B? = -h 2 (t)5t, (27) 
where the index a runs over the SU(2) internal space, while the index i denotes the three spatial directions in the 



spacetime manifold. The use of a cosmic triad parametrized as in Eqs. (271, allows them to include both fields, 
electric and magnetic, and therefore a novanishing invariant G, then an energy density of the form of Eq. ^ with 
X = 1 + 3a" 4 can be derived (cf. Eqs. (17) and (38) in [2D]). 

3) In [6] the isotropization of F and G is done by averaging over a spatial volume, when considered as a classical 
contribution of the radiation, obtaining an effective equation of state for the Born-Infeld field, given by Eq. ([5| , with 
X = 1 + F/(2f3 2 ), for the vanishing electric field case. In the conventions of [BJ, it corresponds to at = (vanishing 
electric field) and 1 = (G — 0), and identifying the coupling constant as A/2 — ► /3 2 , cf. Eqs. (10) in [BJ. 

4) Finally the effective energy density given by Eq. ^ includes the study by Labun and Rafelski |J] ; they consider 
the classical nonlinear Born-Infeld lagrangian as an effective potential, V e &(S, P), that is a nonlinear function of the 
electromagnetic invariants. They obtained the BI energy density given, in our notation, by 

> , j ^»tf- t «V Ft,? ' l '- fl ' t " B ' t, '» (28) 
4 ivX 

cf. Eq. (16) in [9]; that for a magnetic universe E = 0, becomes Eq. ([H]), with X = 1 + F/{2[5 2 ). With the following 
identifications, 

V cS ->L; 4S->F; P -> G; M A -» /3 2 

Our cosmological tests can be applied to their treatment, in what Born-Infeld is concerned. The estimation derived 
from our analysis for the scaling parameter M 4 introduced in [pj is M = 1.127 x 10 _3 eV. 



V. CONCLUSIONS 



Nonlinear electromagnetic fields are very reasonable to assume for inflationary epochs (high energies, anisotropy), 
however in late universes it is not well understood what its origin could be, the conformal symmetry breaking being 
an appealing NLED feature. 

In this paper we have used SNIa, Hubble parameter and GRBs data to probe a cosmological model consisting of 
a homogeneous and isotropic geometry coupled to an electromagnetic field governed by a nonlinear lagrangian of the 
Born-Infeld type. 

The most recent compilation of Supernovae, Union2.1, and the update of the Hubble parameter dataset, provide 
us an improvement in the analysis of cosmological constraints which is reflected in a good statistics allowing us more 
confidence in these results than those obtained from one of the most reliable GRBs dataset. From Supernovae and 
Hubble parameter we found that the BI matter abundance is of Qbi = 0.037 that turns out to be not very significative 
in comparison with the contribution of the cosmological constant which clearly drives the present cosmic acceleration. 

On the other hand, our results from the analysis for the parameter of state w(z) shows a preferrence for quintessence, 
as w(z) > — 1. This result is even clearer when the test is done with GRBs. Moreover the abundance for BI matter 
turns out to be greater for far epochs with redshifts in ranges z ~ 6 — 8, although this sample has a greater dispersion 
in the data in comparation with SNela, and consequently the % 2 = 2.254, points to use of GRBs data always combined 
with another independent probes. 

Moreover, from the found abundance of BI field ftsi = 0.037, we get an estimate of the present Born-Infeld field 
f3 = 1.17 x 10~ 2 Volt/cm, as electric field or f3 = 3.9 x 10 _5 Gauss, that is like ten times the typical strength of magnetic 
fields in galaxies. The corresponding present BI energy density is of the order of (3 2 — 1.618 x 10~ 12 cV 4 . As for the 
BI component the energy density we obtained is four hundredth the critical density, pbi = 0.04p cr i t ~ 10~ 31 gr/cm 3 , 
or pbi = 2.25 x 10 2 eV/cm 3 , that is hundred times the energy density of the CMB radiation; while the estimation for 
the BI coupling constant is g = 0.06Hq in units 1/sec 2 , g = 3.2 x 10 -35 1/sec 2 . 
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To discard dark energy in favor of a vector or tensor field, it should not only accelerate the Universe, but do it at 
the right time and by the right amount, according to observational data. As a result of our tests of the Born-Infeld 
field, we can discard this kind of NLEM field as producing the dark energy effect. However, from Eqs. pi), it can be 
found that 

p=-p- % -(E 2 + B 2 )L F , (29) 
so the possibility exists of an electromagnetic field to account for dark energy or at least to alleviate that problem. 
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